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ABSTRACT 

Swift J1822.3— 1606 was discovered on 2011 July 14 by the Swift Burst Alert Telescope following 
the detection of several bursts. The source was found to have a period of 8.4377 s and was identified 
as a magnetar. Here we present a phase-connected timing analysis and the evolution of the flux and 
spectral properties using RXTE, Swift, and Chandra observations. We measure a spin frequency of 
0.1185154343(8) s"^ and a frequency derivative of -4.3 ± 0.3 x 10"^^ at MJD 55761.0, in a timing 
analysis that include significant non-zero second and third frequency derivatives that we attribute 
to timing noise. This corresponds to an estimated spin-down inferred dipole magnetic field of _B ~ 
5 X 10^^ G, consistent with previous estimates though still possibly affected by unmodelled noise. 
We find that the post-outburst 1-10 keV flux evolution can be characterized by a double-exponential 
decay with decay timescales of 15.5 ± 0.5 and 177 ± 14 days. We also fit the light curve with a crustal 
cooling model which suggests that the cooling results from heat injection into the outer crust. We find 
that the hardness-flux correlation observed in magnetar outbursts also characterizes the outburst of 
Swift J1822.3— 1606. We compare the properties of Swift J1822.3— 1606 with those of other magnetars 
and their outbursts. 

Subject headings: pulsars: individual (Swift J1822.3— 1606) — stars: neutron — X-rays: bursts — 
X-rays: general 



I 1. INTRODUCTION 

Over the past two decades, sever al new classe s of neu- 
tron stars have been discovered fsee lKaspill2010l for a re- 
view). Perhaps the most exotic is that of the magnetars, 
which exhibit some highly unusual properties, often in- 
cluding violent outbursts and high persistent X-ray lumi- 
nos ities that exceed their spin - down powers (for reviews 
• see I Woods fc ThompsonI 120061 : lMereghettill2008| ). These 
' objects, while previously classified as anomalous X-ray 
pulsars (AXPs) and soft gamma repeaters (SGRs), are 
now generally accepted as a unified class of neutron stars 
powered by the decay o f ultra-strong magnetic fields (e.g. 
; [Thompson eTaIll200l . 

To date, there are roughly two dozen magnetars and 
candidates observedQ with spin periods between 2 and 
' 12 s, and high spin-down rates that generally suggest 
- dipole g-fie lds of order 10 ^'^ to 10^^ G (except SGR 
' 0418-H5729: lRea et al.ll2010D . Thanks to the Swift satel- 
\ lite, several new magnetars have be en discovered in 
recent years via t heir outbursts (e.g. iRea et al.l 120091 : 
iGogiis et al.ll2010t UTargaltsev et alll2012[ ). Once a new 
source has been identified, long-term monitoring is cru- 
cial to measure its timing properties, and hence to con- 
strain the dipole magnetic field strength. Also, the flux 
evolution following an outburst could provide insights 
into many physical properties, such as the location of en- 
ergy deposition durin g an outburst, the crust thickness 
and heat capacity (see lPons fc Rea''2012l Gumming et al. 
in prep.), or the physics of a hig hl y activ e magnetosphere 
(|Beloborodovl[2009l : Farfrev et'aLllMl) . 

One of the latest additions to the list of magnetars 
is Swift J1822.3— 1606. This source was first detected 

^ pscholzOphysics.mcgill.ca 

^ See the magnetar catalog at 

|http : //www. physics .mcgill . ca/$\sim$pulsar /magnet ar/maln. html 



by Swtft Burst Alert Telescope (BAT) on 2011 July 14 
(MJD 55756) via its bursting activities (jCummings et al] 
[20T1 . It was soon identified as a new m agnetar upon 
the d etection of a pulse period P=8. 4377 s (jGogiis et al.l 
120 lit ). No optical counterpart was found, with 3a U mit 
down to a z-band mag nitude of 22.2 (|Rea et al.ll2011| ). In 
iLivingstone et al.l (|2Qil, hereafter Paper I), we reported 
initial timing and spectroscopic results using follow-up 
X-ray observations from Swift, Rossi X-ray Timing Ex- 
plorer (RXTE), and Chandra X-ray Observatory. We 
found a spin-down rate of P = 2.54 x 10^^^ which im- 
plies a surface dipole magnetic fielc0 B = 4.7 x 10^'^ G, 
the second lowest i?-ficld among magnetars. Using an 
additional 6 mo nths of Swift and XMM-Newton data, 
IRea et al] (|2012f ) present a timing solution and spectral 

analysis. They find a spin-down rate of P = 8.3 x 10"^'* 
which implies a magnetic field of i? = 2.7 x 10^"^, slightly 
lower than that found in Paper I. 

In this paper, we present an updated timing solution, 
and the latest flux evolution using new observations from 
the same X-ray instruments as in Paper I. The addi- 
tional two Chandra and 18 Swift observations provide 
a timing baseline that is over four times longer and al- 
lows a detailed study of the flux decay. We also report 
on an archival ROSAT observation to constrain the pre- 
outburst flux. We discuss the effects of timing noise on 
our timing solution and the properties and implications 
of this outburst within the magnetar model. 

2. OBSERVATIONS 

2.1. Swift Observations 

The Swift X-Ray Telescope (XRT) consists of a Wolter- 
I telescope and an XMM-Newton EPIC-MOS CCD de- 

^ The surface dipolar component of the B-field can be estimated 
hy B = 3.2 X 10i*'(PP)i/2 q 
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Table 1 

Summary of observations of Swift J1822.3— 1606. 



ObsID Mode Obs Date MJD Exposure Days since trigger 

(TDB) (ks) 



Chandra 



12612 


ASIS-S CC 


2011-07-27 


55769.2 


15.1 


12.6 


13511 


HRC-I 


2011-07-28 


55770.8 


1.2 


14.2 


12613 


ASIS-S CC 


2011-08-04 


55777.1 


13.5 


20.5 


12614 


ASIS-S CC 


2011-09-18 


55822.7 


10.1 


66.1 


12615 


ASIS-S CC 


2011-11-02 


55867.1 


16.3 


110.5 


14330 


ASIS-S CC 


2012-04-19 


56036.9 


20.0 


280.4 


ROSAT 


rp500311n00 




1993-09-12 


49242 


6.7 




Swift 


0(1032(1330(11 


PC 


201 1-07-1 5 


55757.7 


1.6 


1.2 




WT 




^^7^8 7 
OO 1 oo. t 


2.0 


2.1 




WT 


901 1-07-1 7 


55759 7 


2.0 


3.1 


00032033005 


WT 


201 1-07-1 9 


55761.1 


0.5 


4.6 


00032033006 


WT 


2011-07-20 


55762 


1.8 


5.5 


00032033007 


WT 


201 1-07-21 


55763 2 


1.6 


6.7 


00032033008 


WT 


201 1-07-23 


55765 8 


2.2 


9.2 


00032033009 


WT 


2011-07-24 


55766 2 


1.7 


9.7 


0003203301 


WT 


201 1-07-27 


55769 5 


2.1 


12.9 


0003203301 1 


WT 


201 1-07-28 


55770 3 


2.1 


13.8 


00032033012 


WT 


201 1-07-29 


55771.2 


2.1 


14.7 


0003203301 3 


WT 


201 1-07-30 


55772 3 


2.1 


15.7 


00032051 001 


WT 


201 1-08-05 


55778 


1.7 


21.5 


00032051 002 


WT 


201 1-08-06 


55779 


1.7 


22.5 


00032051003 


WT 


201 1-08-07 


55780 4 


2.3 


23.9 


00032051 004 


WT 


201 1-08-08 


55781.4 


2.3 


24.8 


00032051 005 


WT 


201 1-08-1 3 


55786 4 


2.2 


29.8 


00032051006 


WT 


2011-08-14 


55787 6 


2.2 


31.0 


00032051 007 


WT 


201 1-08-1 5 


55788 1 

'J'J 1 OO . -L 


2.3 


31.6 


00032051 008 


WT 


201 1-08-1 6 


55789 5 


2.2 


32.9 


00032051009 


WT 


2011-08-17 


55790 3 


2.2 


33.8 


00032033015 


WT 


2011-08-27 


55800.8 


2.9 


44.2 


00032033016 


WT 


2011-09-03 


55807.2 


2.4 


50.6 


00032033017 


PC 


2011-09-18 


55822.7 


4.9 


66.2 


00032033018 


WT 


2011-09-20 


55824.5 


1.5 


68.0 


00032033019 


WT 


2011-09-25 


55829.1 


2.3 


72.6 


00032033020 


WT 


2011-10-01 


55835.1 


2.6 


78.5 


00032033021 


WT 


2011-10-07 


55841.7 


4.2 


85.2 


00032033022 


WT 


2011-10-15 


55849.2 


3.4 


92.7 


00032033023 


WT 


2011-10-22 


55856.2 


2.2 


99.7 


00032033024 


PC 


2011-10-28 


55862.2 


10.2 


105.6 


00032033025 


PC 


2012-02-19 


55976.4 


6.2 


219.8 


00032033026 


WT 


2012-02-20 


55977.0 


10.2 


220.5 


00032033027 


PC 


2012-02-21 


55978.1 


11.0 


221.6 


00032033028 


WT 


2012-02-24 


55981.9 


6.7 


225.4 


00032033029 


WT 


2012-02-25 


55982.8 


7.0 


226.3 


00032033030 


WT 


2012-02-28 


55985.0 


7.0 


228.5 


00032033031 


WT 


2012-03-05 


55991.1 


6.8 


234.5 


00032033032 


WT 


2012-04-14 


56031.1 


4.3 


274.6 


00032033033 


WT 


2012-05-05 


56052.6 


5.1 


296.0 


00032033034 


WT 


2012-05-26 


56073.0 


4.9 


316.5 


00032033035 


WT 


2012-06-17 


56095.5 


5.6 


338.9 


00032033036 


WT 


2012-06-26 


56104.1 


6.2 


347.6 


00032033037 


WT 


2012-07-06 


56114.2 


6.8 


357.6 


00032033039 


WT 


2012-08-17 


56156.1 


4.9 


399.6 


00032033040 


WT 


2012-08-22 


56161.5 


5.0 


405.0 



tector ("B urrows et al.l 12005'). Swift is optimized to pro- 
vide rapid follow-up of gamma-ray bursts and other X- 
ray transients. Following the 2011 July 14 outburst of 
Swift J1822.3-1606, the XRT was used to obtain 46 ob- 
servations for a total exposure time of 175 ks. Data were 
collected in two different modes, Photon Counting (PC) 
and Windowed Timing (WT) . While the former gives full 
imaging capability with a time resolution of 2.5 s, the lat- 
ter forgoes imaging to provide 1.76-ms time resolution by 
reading out events in a collapsed one-dimensional strip. 
For each observation, the unfiltered Level 1 data were 



downloaded from the Swift quicklook archive^. For 
a summary of observations used, see Table [T] The 
standard XRT data reduction script, xrtpipeline, was 
then run using HEASOFT 6.11 and the Swift 20110725 
CALDB. We reduced the events to the barycenter using 
the position of RA = 18*^ 22™ 18", Dec= -16° 04' 26'.'8 
(iPagani et al.ll20lil ). Source and background events were 
extracted using the following regions: for WT mode, a 
40-pixel long strip centered on the source was used to 
extract the source events and a strip of the same size 

* http:/ /swift. gsfc.nasa.gov/cgi-bin/sdc/ql 
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Table 2 

Summary of RXTE observations of Swift J1822.3-1606. 



ObsID Obs Date MJD Exposure Days since 

(TDB) (ks) trigger 



D96048- 


-02- 


-01- 


-00 


2011- 


-07- 


-16 


55758, 


.49 


6, 


.5 


1.96 


D96048- 


-02- 


-01- 


-05 


2011- 


-07- 


-18 


55760, 


,81 


1, 


.7 


4.28 


D96048- 


-02- 


-01- 


-01 


2011- 


-07- 


-19 


55761, 


.57 


5, 


.1 


5.04 


D96048- 


-02- 


-01- 


-02 


2011- 


-07- 


-20 


55762, 


.48 


4, 


.9 


5.95 


D96048- 


-02- 


-01- 


-04 


2011- 


-07- 


-21 


55763, 


.42 


3, 


.3 


6.89 


D96048- 


-02- 


-01- 


-03 


2011- 


-07- 


■21 


55763, 


.64 


6, 


.0 


7.11 


D96048- 


-02- 


-02- 


-00 


2011- 


-07- 


-22 


55764, 


.62 


6, 


.1 


8.09 


D96048- 


-02- 


-02- 


-01 


2011- 


-07- 


-23 


55765, 


.47 


6, 


.8 


9.94 


D96048- 


-02- 


-02- 


-02 


2011- 


-07- 


-25 


55767, 


.60 


3, 


.0 


11.07 


D96048- 


-02- 


-03- 


-00 


2011- 


-07- 


-29 


55771, 


.35 


6, 


.8 


14.82 


D96048- 


-02- 


-03- 


-01 


2011- 


-08- 


-01 


55774, 


.35 


6, 


.9 


17.82 


D96048- 


-02- 


-03- 


-02 


2011- 


-08- 


-04 


55777, 


.85 


1, 


.9 


21.32 


D96048- 


-02- 


-03- 


-04 


2011- 


-08- 


-04 


55777, 


.92 


1, 


.8 


21.39 


D96048- 


-02- 


-04- 


-00 


2011- 


-08- 


-07 


55780, 


.49 


6, 


.9 


23.96 


D96048- 


-02- 


-04- 


-01 


2011- 


-08- 


-09 


55782, 


.58 


6, 


.5 


26.05 


D96048- 


-02- 


-04- 


-02 


2011- 


-08- 


-11 


55784, 


.97 


3, 


.7 


28.44 


D96048- 


-02- 


-05- 


-02 


2011- 


-08- 


-12 


55785, 


.03 


3, 


.3 


28.50 


D96048- 


-02- 


-05- 


-00 


2011- 


-08- 


-15 


55788, 


.05 


5, 


.9 


31.52 


D96048- 


-02- 


-05- 


-01 


2011- 


-08- 


-16 


55789, 


.96 


6, 


.0 


33.43 


D96048- 


-02- 


-06- 


-00 


2011- 


-08- 


-21 


55794, 


.46 


6, 


.6 


37.93 


D96048- 


-02- 


-07- 


-00 


2011- 


-08- 


-26 


55799, 


.61 


6, 


.8 


43.1 


D96048- 


-02- 


-08- 


-00 


2011- 


-09- 


-06 


55810, 


.38 


6, 


.0 


53.8 


D96048- 


-02- 


-10- 


-00 


2011- 


-09- 


-16 


55820, 


.24 


6, 


.7 


63.7 


D96048- 


-02- 


-10- 


-01 


2011- 


-09- 


-22 


55826, 


.18 


5, 


.6 


69.6 


D96048- 


-02- 


-09- 


-00 


2011- 


-09- 


-25 


55829, 


.38 


6, 


.2 


72.8 


D96048- 


-02- 


-11- 


-00 


2011- 


-10- 


-01 


55835, 


.90 


7, 


.1 


79.4 


D96048- 


-02- 


-12- 


-00 


2011- 


-10- 


-08 


55842, 


.23 


5, 


.9 


85.7 


D96048- 


-02- 


-13- 


-00 


2011- 


-10- 


-15 


55849, 


.67 


5, 


.6 


93.1 


D96048- 


-02- 


-14- 


-00 


2011- 


-10- 


-29 


55863, 


.11 


6, 


.7 


106.6 


D96048- 


-02- 


-16- 


-00 


2011- 


-11- 


-13 


55878, 


.90 


5, 


.9 


122.4 


D96048- 


-02- 


-17- 


-00 


2011- 


-11- 


-20 


55885, 


.21 


6, 


.0 


128.7 


D96048- 


-02- 


-15- 


-00 


2011- 


-11- 


-28 


55893, 


.18 


6, 


.7 


136.6 



positioned away from the source was used to extract the 
background events. For PC mode, a circular region with 
radius 20 pixels was used for the source region and an 
annulus with inner radius 40 pixels and outer radius 60 
pixels was used as the background region. For the first 
PC mode observation (00032033001), a circular region 
with radius 6 pixels was excluded to avoid pileup. For 
the subsequent PC mode observation (00032033017), a 
region with radius 2 pixels was excluded. We estimate 
the maximum pileup fraction of the remaining PC obser- 
vations be less than 5%. 

For the spectral analysis, Swift ancillary response files, 
which provide the effective area as a function of energy, 
were created using the FTOOL xrtmkarf and the spec- 
tral redistribution matrices from the Swift CALDB were 
used. 

2.2. RXTE observations 

The RXTE Proportional Counter Array (PCA) com- 
prised five proportional counting units and provided 
a large collecting area and high timing precision 
(jJahoda et al.1 l2006l ). We downloaded 32 observations 
from the HEASARC archive spanning an MJD range 
from 55758 to 55893, for a total of 174 ks of integra- 
tion time. The data were collected in GoodXenon mode 
which records each event with l-/is time resolution. The 
observations are summarized in Table [2j 

We selected events in the 2-10 keV energy range (PCA 
channels 6-14) from the top xenon layer of each PCU for 
our analysis, to maximize signal-to-noise ratio. The data 
from all the active PCUs were then merged. If more 



than one observation occurred in a 24-hr period, the ob- 
servations were combined into a single data set. Photon 
arrival times were adjusted to the solar system barycen- 
ter using the same position as the for Swift data. Events 
were then binned into time series with resolution 1 /32 s 
for use in the following analysis. 

2.3. Chandra observations 

Following the outburst, we triggered our ToO program 
with the Chandra X-ray Observatory. The telescope on- 
board Chandra has an effective area ~3 times larger than 
that of Swift XRT, when used with the ACIS detector 
in continuous clocking (CC) mode. This mode has a 
time resolution of 2.85 ms and sensitivity between 0.3 
and lOkeV0- Five ACIS CC-mode observations were ob- 
tained between days 13 and 281 after the outburst, with 
exposures ranging from 10 to 20 ks. The observation pa- 
rameters are summarized in Table [TJ For imaging pur- 
poses, we also processed a short (1.2 ks) archival Chandra 
HRC-I observation taken 14 days after the outburst. 

All Chandra data were processed using CIAO 4.3 with 
CALDB 4.4.6. We extracted the source events with a 
6"-long strip region, and the remainder of the collapsed 
strip (^7' long), excluding the region within 1' of the 
source in order to minimize any contamination from the 
wings of the PSF, was used for the background. We 
restricted the timing analysis to events between 0.3 and 
8keV. Photon arrival times were corrected to the solar 
system barycenter. The source spectrum was extracted 
using the tool specextract. 

2.4. ROSAT observation 

The only existing X-ray image that covers the field 
prior to the out burst is a 6.5-ks ROSAT PSPC 
(j AschenbachI 1 1 985l ) observation of the nearby Hii region 
M17 (Omega Nebula, G15.1-0.8). The observation has 
a time resolution of 130 ms. We downloaded the filtered 
event list from the HEASARC data archiv^ and carried 
out the analysis using FTOOLS. 

3. ANALYSIS & RESULTS 
3.1. Imaging 

Figure [1] shows the ROSAT and Chandra images. 
Swift J1822.3— 1606 is the only source detected in the 
Chandra HRC image and its radial profile is fully con- 
sistent with that of a model PSF. Hence, there is no 
evidence for any surrounding nebula or dust scattering 
halo. We find a source position of RA= 18*^ 22™ 18.06", 
Dec= —16° 04' 25'.'55 from the HRC image which iscon- 
sistent with the XRT position from lPagani et aP ()2011[ ) 
used above. We assume an error radius of 0.6" which 
is the uncertainty in the absolute astrometry of Chan- 
dra for a 90% confidence interval Q- In the ROSAT im- 
age, an unresolved source is clearly detec ted at the posi- 
tion of the magnetar, as first reported bv lEsposito et all 
(jUIlli). Since the Chandra image shows no other bright 
X-ray sources in the field, we take this source to be 
Swift J1822. 3-1606 in quiescence. Using a 4' x 2' el- 
liptical aperture, we obtain 113 ± 11 total counts in 0.1- 

^ |http : 1 1 cue . h arvard ■ edu/proposer/P OG/html/] 

^ Ihttp : //heasarc . gsf c . nasa . gov/W3Bro« se/ 1 

^ According to http: / /cxc .harvard. edu /cal/ASPECT/celmon/] 
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2.4 keV range, of which 48 ± 7 counts are due to back- 
ground. Finally, we note that the diffuse X-ray emission 
~20' southwest of the magnetar is from M17, which con- 
tains the young stella r cluster NGC 6618 with over 100 
OB stars ' (|Lada et al.lll991| ). 



Figure 1. Left: ROSAT image of the field of Swift J1822. 3-1606 in 0.1-2.4keV range. The position of Swift J1822. 3-1606 is marked 
by the cross, and the lines indicate the Chandra HRC observation field of view. The large-scale diffuse emission is the Galactic Hll region 
M17. Right: 1.2ks Chandra HRC exposure of Swift J1822. 3-1606. Swift J1822.3-1606 is the only source detected. 

each photon folded at the best ephemeris of the pulsar. 
The likelihood distribution that results then describes 
the probability density for the average pulse arrival time. 
A TOA can be calculated from the optimal phase off- 
set. We estimated TOA uncertainties by simulating one- 
hundred sets of events drawn from the pulse profile of the 
observation and measured an offset for each set using the 
ML method. The standard deviation of the simulated off- 
set distribution was then taken as the TOA uncertainty. 
The ML derived TOAs were consistent with those de- 
rived for Paper I using the cross-correlation method. 

Timing solutions were then fit to the TOAs using 
TEMPC^. Three solutions, one with a single freqency 
derivative, one with two frequency derivatives and one 
with three frequency derivatives, are given in Table [3l 
The top panel of Figure [2] shows the timing residuals 
with just V and v fitted (Solution 1), the middle panel 
shows the residuals with i) also fitted (Solution 2), and 
the bottom panel shows the residuals with i) and v also 
fitted (Solution 3). Solution 1 is a poor fit with a xl/v 
of 5.02/72. This is likely due to timing noise, a com- 
mon phenomenon in young neutron stars including mag- 
netars (e.g. Dib et al. 2008; Livingstone & Kaspi 2011). 
The best-fit v and v values for Solution 1 imply a sur- 
face dipolar magnetic field of B = 2.43 ± 0.03 x 10^^ G. 
Solution 2, with a significant non-zero v, gives a bet- 
ter fit with a y^/v of 1.94/71. An i^-test gives a 
probability of 2 x 10^^^ that the addition of a second 
derivative does not significantly improve the fit. The 
surface dipolar magnetic field implied by Solution 2 is 
B = 3.84 ± 0.08 X 10^3 G. Solution 3, with a significant 
non-zero v , provides still a better fit than Solution 2 
with a x^/i^ of 1.44/70. An F-test gives a probability 
of 3 X 10^^ that the addition of a third derivative does 
not significantly improve the fit. The best-fit parameters 
from Solution 3 imply a surface dipolar magnetic field of 
B = 5.1 ±0.2 X 10" G. 

Note that the fit is heavily influenced by the very high 
quality Chandra TOAs. However, omitting them and 



3.2. Timing Analysis 

3.2.1. Spm Evolution 

Barycentered events were used to derive a pulse time- 
of-arrival (TOA) for each Swift WT mode, RXTE, and 
Chandra observation. For the RXTE observations, 
events were binned into time series with 31.25-ms res- 
olution. The time series were then folded with 128 phase 
bins using the ephemeris from Paper L A TOA was then 
measured from each profile by cross-correlation with a 
template profile. We verified that the RXTE pulse pro- 
files were consistent with each other except in one iso- 
lated observation which was handled accordingly (see 

Section 13X11) • 

For Swift and Chandra observations, TOAs were ex- 
tracted using a Maximum Likelihood (ML) method, as 
it yields more accurate T OAs than the traditional cross- 
correlation technique (see [Livingstone et al.ll2009l ). This 
method was not used for the RXTE observations as their 
high number of counts (due to the large collecting area 
and background count rates of the PCA) make the ML 
method computationally expensive. The ML method for 
measuring TOAs requires a continuous model of the tem- 
plate pulse profile for which we used a Fourier model. 
The discrete Fourier Transform of the binned template 
profile was first calculated. The template was then fitted 
by /(^) = Sj=o '^J^*^'^'''^ where aj is the Fourier coef- 
ficient for the j*^ harmonic, and (p the phase between 
and 1. The number of harmonics used was optimized to 
account for the features of the light curve while ignoring 
small fluctuations caused by the flnite number of counts. 
For Swift J1822.3— 1606, we used flve harmonics to derive 
the TOAs. 

For each observation, a probability or likelihood for 
a grid of trial offsets, (/)off, can be calculated using 
P{.4>os) = 11^=0 f{.4'i - '/'off), where 0^ is the phase of 



* http:/ /www. atnf.csiro.au/people/pulsar/tempo/ 
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including only TOAs from Swift and RXTE still yields 
significant second and third derivatives and an implied 
B-field oiB = 4.8 ± 0.2 x 10" G which is consistent with 
that of Solution 3. The above-quoted uncertainties in 
B and other derived quantities in Table [3] reflect only 
the statistical uncertainties in u and its derivatives and 
do not include any contributions from the simplified as- 
sumptions in the standard formulae used to determine 
such quantities. Note that even with the addition of 
highly significant second and third derivatives, Solution 
3 still does not provide an adequate fit. Adding addi- 
tional derivatives reduces the with marginal signifi- 
cance and results in larger values of the spin-down rate 
and hence B. For example, including a fourth frequency 
derivative does not result in significant improvement in 
{xl/v = 1.31/69) and yields B = 6.0 x lO^^ G. 
To search for pulsations in the ROSAT observation, 
we applied a barycenter correc tion to the event arr ival 
times, then used the Z"^ test (jBuccheri et al.l IT983D to 
search for pulsations. We searched in the frequency range 
from zero to 3.8 kHz in steps of 1.3 /iHz, oversampling the 
independent Fourier spacing by a factor of 10; however, 
we found no significant signal. By simulating a pulsar 
with a background subtracted count rate of that of the 
ROSAT observation, we find that the pulsar would be 
undetectable even with a pulsed fraction of 100%, there- 
fore we cannot constrain the pulsed fraction. 

3.2.2. Pulse Profile Analysis 

Here we search for time and energy variability in the 
pulse profile of Swift J1822.3-1606 using the RXTE, 
Swift, and Chandra observations. We created pulse pro- 
files for each RXTE observation for energy ranges of 2 - 
6keV, 6 - lOkeV (with photons selected from only the 
top xenon layer), 10 - 15keV, 15 - 20keV, 20 - 4dkeV, 
and 20 - 60keV (with photons selected from all three 
xenon layers) using the Solution 2 ephemeris. For the 
Chandra data we produced pulse profiles with the en- 
ergy ranges of 0.5 - 6keV, 0.5 - 2keV and 2 - 6keV. 
For the Swift data we created 0.5 - lOkeV profiles, using 
only WT mode observations as PC mode does not have 
sufficient time resolution. 

As in Paper I, we searched for time variability in pulse 
profiles but found that all the RXTE profiles are consis- 
tent with the template in each case except for the one 
profile from the very first observation after the outburst 
(Obsid D96048-02-01-00). The difference is due primarily 
to the off-pulse feature which had slightly different struc- 
ture between the template profile and the first RXTE 
observation. We therefore did not use D96048-02-01-00 
in the timing analysis. 

The Chandra profiles, however, do show evidence 
for low-level variability. Figure [3] shows the 0.5 - 
6 keV pulse profile for each Chandra observation of 
Swift J1822.3— 1606. We produced residuals between 
each pair of Chandra profiles by normalizing each profile 
and taking the difference between each normalized pair. 
A comparison of the profile residuals between each set 
of profiles shows that there is significant low-level evo- 
lution of the small 'pulse' that precedes the main pulse. 
The main pulse does not exhibit any significant variation. 
The most significant variability is that between the first 
(MJD 55769) and last (MJD 56036) Chandra observa- 
tion. The residuals between those two profiles have a 



Table 3 

Spin Parameters for Swift J1822. 3-1606. 



Parameter 


Value 


Dates (Modified Julian Day) 


55759 - 56161 


Epoch (Modified Julian Day) 


55761.0 


Number of TOAs - RXTE 


31 


Number of TOAs - Swift 


40 


Number of TOAs - Chandra 


5 


Solution 1 - one frequency derivative 


u (s-l) 


0.1185154253(3) 


U (s-2) 


-9.6(3) X lO-i*^ 


RMS residuals (ms) 


52.2 




5.02/72 


B{G) 


2.43(3) X lO" 


E (erg s -"-) 


4.5(1) X 10^'' 


Tc (kyr) 


1963(51) 



Solution 2 - two frequency derivatives 



V (s-l) 


0.1185154306(5) 


!> (s-2) 


-2.4(1) X 10-1^ 


i> (s-3) 


1.12(8) X 10-22 


RMS residuals (ms) 


32.2 




1.94/71 


B{G) 


3.84(8) X 10" 


E (erg s 1) 


1.12(5) X 10^1 


Tc (kyr) 


784(33) 



Solution 3 - three frequency derivatives 



u (s-l) 


0.1185154343(8) 


U (s-2) 


-4.3(3) X 10-1^ 


V (s-3) 


4.4(6) X 10-22 


■v (s-4) 


-2.2(4) X 10-29 


RMS residuals (ms) 


27.5 




1.44/70 


B{G) 


5.1(2) X 10" 


E (erg s 1) 


2.0(2) X 10^1 


Tc (kyr) 


442(33) 



Note. — E 



rrors arc 



formal Icr TEMPO uncertainties. 



of 16.8 for 63 degrees of freedom. We note however that 
these low-level variations in the smaller component likely 
do not have a significant impact on the timing analysis, 
since the TOA extraction is heavily weighted toward the 
unchanging primary component. Indeed our simulations 
of the effects of such low-level profile variations on the 
TOAs (see below) strongly support this conclusion. 

The Swift profiles also show evidence for low-level vari- 
ability. As above, we produced residuals between each 
pair of profiles and calculated a for the null hypoth- 
esis. The measured values are not consistent with a x^ 
distribution, so there is significant variation between pro- 
file pairs. A closer look at the residuals shows that the 
variation is due primarily to the small interpulse, as in 
the Chandra data. 

To investigate the dependence of pulse morphology on 
energy, we created a single high significance profile by 
aligning and summing individual profiles for each energy 
range. Figure |4] shows a summary of the results, with 
the summed profiles for 0.5 - 2keV (top panel, with 64 
phase bins, Chandra data), 2 - 6keV (middle panel, with 
64 phase bins, Chandra data) ,6-10 keV and 10-15 keV 
(bottom two panels, RXTE data, with 64 and 16 phase 
bins, respectively). No pulsations were detected above 
15keV with the PCA. We then calculated residuals be- 
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55750 55800 55850 55900 55950 56000 56050 56100 56150 56200 

MJD 

Figure 2. Timing residuals of Swift J1822.3— 1606. The top panel shows the residuals for the timing solution with one frequency derivative 
(Solution 1); see Table |3] The middle panel shows the residuals with the addition of a second derivative (Solution 2). The bottom panel 
shows the residuals with three frequency derivatives fitted. Swift WT mode data are represented by crosses, red circles denote RXTE 
observations and Chandra observations are shown as blue triangles. 



3 




Figure 3. Pulse profiles with 64 bins from each of the four Chan- 
dra observations. The profiles are for the energy band 0.5-6 keV. 

tween pairs of profiles, and calculated xt values of the 
resulting residuals in order to identify energy dependence 
of the pulse morphology. The most significant variability 




Phase (2 cycles) 

Figure 4. Pulse profiles for Swift J1822.3— 1606 from Chandra 
and RXTE data for four energy ranges: 0.5-2 keV and 2-6 keV 
{Chandra data, 64 phase bins), 6-10 kcV and 10-15 keV (RXTE 
data, 64 and 16 phase bins). 



is between the 0.5 - 2keV Chandra profile and the 6 - 
lOkeV RXTE profile. This can been seen in Figure |3] as 
a change in the phase of the interpulse, arriving later for 
higher energies. For this profile pair, the xt of the resid- 
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uals is 46.2 (for 28 degrees of freedom), excluding the 
null hypothesis. The interpulse variability causes signifi- 
cant differences between each pair of profiles, except the 
10 - 15keV profile, likely because of the lower statistics 
of the latter. 

3.2.3. Comparison to Previously Reported Results 

iRea et all (|2012[ ) present a timing solution with a spin- 
down implied magnetic field of B — 2.7 x 10^'^G. Their 
data set is similar to ours, although they use proprietary 
XMM-Newton data whereas we use proprietary Chandra 
data and our data set includes seven additional Swift 
observations. Their timing solution is similar to our So- 
lution 1. They, however, do not find a significant second 
frequency derivative. A possible cause of this discrep- 
ancy could be the difference in TOA extrac tion meth- 
ods. I nstead of using a pulse profile template, IRea et ahl 
()2012t ) fit the folded profile for each observation with two 
sine functions with periods equal to the fundamental and 
the first harmonic of the pulse period. They then assign 
the ascending node of the fundamental sine function as 
the time-of-arrival of the pulse. This method was used 
to attempt to account for pulse-profile changes. We im- 
plemented this method and derived an additional set of 
TOAs to compare to our ML derived TOAs. We found 
that the sine-model derived TOAs provided similar tim- 
ing solutions as our ML TOAs and the addition of a 
second frequency derivative did significantly improve the 
fit, reducing the xl/i' from 7.91/72 to 2.72/71. The ad- 
dition of a third derivative in this case, results in only 
marginal improvement with a x^/*^ of 2.47/70. If we 
limit our d atase t to the Swift and RXTE data used in 
IRea et al.l (|20T2| ). we find that the addition of a second 
frequency derivative is not necessary, which is consistent 
with their findings. 

In order to investigate the effects of pulse profile 
changes on both TOA extraction methods, we simu- 
lated pulse profiles with an unchanging (other than noise) 
primary component and a varying secondary compo- 
nent, as is observed in the pulse profile evolution of 
Swift J1822.3— 1606. We modelled the profile using two 
gaussians and modified the amplitude of the smaller 
gaussian in order to vary the secondary component. We 
found that for both the sine-model and the ML methods, 
as the amplitude of the secondary component was varied, 
the measured phase offsets varied by less than their un- 
certainties. Hence we conclude that the observed pulse 
profile variations do not have an appreciable affect on the 
TOA determination, independent of which TOA extrac- 
tion method was used. 

3.3. Spectral Analysis & Flux Evolution 

Spectral models were fit to the Swift, Chandra and 
ROSAT data using XSPECfl vl2.7. The quiescent flux 
of Swift J1822.3— 1606 was determined by first extracting 
the source spectrum from the ROSAT data, then fitting 
it with an absorbed blackbody model. The absorption 
column density A'h was fixed during the fit at the best- 
fit value (4.53 x lO^^cm"^) determined from the Swift 
and Chandra spectra (see below). We obtained a qui- 
escent blackbody temperature of kT = 0.12 ± 0.02 keV 

^ http://xspoc.gfsc.nasa.gov 
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Figure 5. The spectral properties of Swift J1822.3— 1606 follow- 
ing its outburst. Grey crosses represent Swift WT mode data, open 
circles are Swift PC mode data, and blue triangles denote Chandra 
observations. The determination of t he blackbody radius assumes 
a distance of 1.6 kpc (see Section |4.3| I. 




100 150 200 250 300 350 
Days from BAT Trigger (MJD 55756.5) 



Figure 6. Flux evolution of Swift J1822. 3-1606 following its 
outburst. The top panel shows the 1— lOkeV absorbed flux as it 
decays. The middle panel shows that 2-10 keV pulsed count rate. 
The RXTE /PCA and Chandra count rates have been scaled to the 
Swift count rates by fitting overall scaling factors in the double- 
exponential fit. The bottom panel shows the pulsed fraction from 
Swift and Chandra observations. In all panels, Swift WT mode 
data are shown as grey crosses, Swift PC mode data are represented 
by open circles, RXTE data are shown as red circles and data 
from Chandra are shown as blue triangles. The dotted lines show 
the result of the double-exponential fit to the total and pulsed 
flux decays, and the dashed line shows the results of the single- 
exponential fit. 

and a radius of SljC^i.gkm, where di.g is the distance to 
the source in units of 1.6 kpc, the estimated distance as 
discussed in Section 14.31 We found an absorbed fiux of 
9tg° X 10-i^ergcm~2s~i in the 0.1-2.4keV range. 
The Chandra and Swift spectra were grouped with 
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Table 4 

Models of the Flux Evolution of Swift J1822. 3-1606. 



Model Ti (days) (days) Fq (erg cm ^s ^) X^/'^ 

Single Exponential 23.8 ± 0.5 - 3 X 10"^" (fixed) 20.1/43 

Single Exponential 19.5 ± 0.4 - 4.7 ± 0.2 X lO^^^ 4,48/42 

Double Exponential 15.5 ± 0.5 177 ± 14 3 X lO'i"! (gxed) 2.17/41 

Double Exponential 9±1 39 ± 3 4.0 ± 0.2 X lO"^^ 1.06/40 



a minimum of 100 and 20 counts per bin, respec- 
tively. The spectra were fitted jointly to a photoelectri- 
cally absorbed blackbody model with an added power- 
law component. The model was fit with a single iVjj 
using the XSPEC phabs model assuming abundances 
fromTAndcrs & Grcvcssc (1989) and photoelectric cross- 
sections from lBalucihska-Church fc McCammonI (|1992i ). 
All the other parameters were allowed to vary from ob- 
servation to observation. 

The blackbody plus power-law model provided an ac- 
ceptable fit to the Swift and Chandra data. The model 
had -y^^lv of 1.07/5451 and a best-fit iVn of 4.53 ±0.08 x 
lO^^cm"^. Figure [5] shows the evolution of the post- 
outburst spectral parameters. The spectrum is softening 
following the outburst as the 1-10 keV absorbed flux 
(top panel of Figure [6]) decays: the photon index of the 
power-law component increases and the temperature of 
the blackbody decreases. This is clear from the high- 
quality Chandra data alone but is also apparent in the 
Swift data, and is consistent with the behavior of other 
magnetars post-outburst (see Section 14. 2p 

Figure [6] shows the flux decay and pulsed fraction evo- 
lution following the Swift J1822.3— 1606 outburst. In Pa- 
per I, we showed that both a double exponential and 
an exponential model provided acceptable fits to both 
the total and pulsed flux decays, whereas a power-law 
decay model was excluded. Here we fitted a double- 
exponential decay and an exponential decay to the to- 
tal and pulsed flux evolutions. The exponential de- 
cay model is described by F{t) — Fp exp^^*"*")/"^ -l-F, 
where t is in MJD, Fp is the peak absorbed flux, 
Fq is the quiescent flux, tg is the time of the BAT 
trigger in MJD, and r is decay timescale in days. 
The double-exponential decay is described by F{t) = 
Fiexp-(*-*o)/^i -fF2exp-(*-*o)/^2 _^Fg where to is the 
time of the BAT trigger in MJD, Fi and F2 are the ab- 
sorbed fluxes at ^0 of each exponential component and ti 
and T2 are the decay timescales in days of each compo- 
nent. For both models we fit the data both by using the 
quiescent fiux as a free parameter, and by using a fixed 
quiescent flux = 3 x 10~^^ ergcm"^ s~^. This is the 
approximate 1-lOkeV flux assuming the 0.1-2.4keV flux 
and spectral model from the ROSAT observation. Table 
m shows the results of these fits. With xl/^ of 1.06/40, 
the double-exponential decay with a free quiescent flux 
provides the best fit to the data. However, the best-fit 
value for the quiescent fiux, 4.0±0.2x 10~^^ergcm~^ s~^, 
is more than two orders of magnitude higher than the 
quiescent flux implied by the ROSAT observation. We 
therefore take the double-exponential decay with the 
fixed quiescent fiux as the best model of the total ab- 
sorbed fiux decay with timescales of ti = 15.5 ± 0.5 
days and T2 — 177 ± 14 days. In both cases, the single- 



exponential fit is much worse than the double exponen- 
tial. 

In Figure [SI the absorbed flux measured with Chan- 
dra is always larger than that with Swift at a similar 
epoch by about 10-15%. The discrepancy could be at- 
tribute d to insufficient cross-c alibration between instru- 
ments. iTsuiimoto et al.l (|2011| ) found that flux measure- 
ments from different X-ray telescopes could differ by as 
much as 20%, and that Chandra appears to give a higher 
flux than others, as well as a harder photon-index, which 
is consistent with our findings. 

To determine the 2-10 keV pulsed count rate from 
Swift, RXTE, and Chandra observations, the barycen- 
tered events were folded using the Solution 2 ephemeris 
in Table [3] with 16 phase bins. For the RXTE obser- 
vations, only data from the first xenon layer of PCU2 
were used. Both PCUO and PCUl lost their propane 
layers and there is minimal PCU3 and PCU4 data for 
this source. The inclusion of only a single detector in the 
analysis reduces instrumental biases. For observations 
from all three telescopes, the pulsed count rate was then 
measured from each folded profi le using a RMS method 
as described in lDib et al.l (|2008[ ). where the pulsed count 
rate, F, is given by: 



F = 



2E[K + ^i)-«+<)] 



(1) 



k=l 



where is the even Fourier component and is equal 
to [l/N] cos{2Trki/N), Ga^^ is the uncertainty in 

afc, 6fe is the odd Fourier component and is equal to 
(1/^) Er=iP^sin(2^fcz/iV), Oht, is the uncertainty in 6fe, 
i is an index over phase bins, iV is the total number of 
phase bins, is the count rate in the zth phase bin, and 
n is the maximum number of Fourier harmonics used. 
In this case, n = 5. This technique is equivalent to the 
simple RMS formula, F = {l/VN)[J2f^^{pi-p)^Y/^, ex- 
cept only statistically significant Fourier components are 
used and the upward statistical bias is removed by sub- 
tracting the variances (Dib et al.' 2008). For Swift and 
Chandra observations, pulsed fractions were determined 
by dividing the pulsed count rate by the total count rate. 

The middle panel of Figure [S] shows 
the 2-10 keV pulsed-flux evolution of 
Swift J1822.3— 1606. The pulsed count rates mea- 
sured by each instrument depend on the different 
instrumental responses. The RXTE PGA and Chandra 
pulsed count rates were therefore scaled to the Swift 
WT mode values by including factors between each 
data set as free parameters in the double and single- 
exponential fits. For the pulsed-flux evolution, the 
double-exponential fit also provided the best fit with 
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xi/y of 5.15/74 and decay timescales of n = 15.3 ± 0.2 
days and T2 = 182 ± 6 days. The exponential model 
had a xi/v of 60.5/76 with a best-fit decay timescale of 
25.1 ± 0.2 days. This is the opposite of what we found 
in Paper I where the exponential model was a better fit 
to the data available at the time. 

3.4. X-ray Bursts 

To search for X-ray bursts in RXTE data of 
Swift J1822.3— 1606, we created a time series for each 
active PCU from GoodXenon data for each observation, 
selecting events in the 2-20 keV range (PCA channels 6- 
24) and from all three detection layers (the same energy 
range as selected for similar searches for X-ray bursts 
from magnetars, e.g. IGogiis et al.l 120011 : IGavriil et al.l 
l2004f ). For the Swift observations, binned time series 
were made for each Good Timing Interval (GTI) in an 
observation. For both Swift and RXTE, time series were 
made at 15.625-ms, 31.25-ms, 62.5-ms and 125-ms time 
resolutions to provide sensitivity to bursts on a hierarchy 
of time scales. 

Bursts were identified by comparing the count rate in 
the i*^ bin t o the average count rate as described in 
IGavriil et al.l (2004). Because the background rate of 
the PCA typically varies over a single observation, we 
calculated a local mean around the i*'* bin for RXTE. 
For Swift data, a mean was calculated for each GTI. We 
then compared the count rate in the i'^ bin to the mean. 
If the count rate in a single bin was larger than the lo- 
cal/GTI average, the probability of such a count rate 
occuring by chance was calculated. For RXTE data, the 
probability of the count rate in the corresponding bin in 
the other active PCUs was also calculated (whether or 
not the count rate in that bin was greater than the local 
average). If a PCU was off during the bin of interest, 
its probability was set to 1. We then found the total 
probability that a burst was observed, by multipying the 
probabilities for each PCU together. If the total proba- 
bility of an event was Pi, tot < O.Ol/TV (where TV is the 
total number of time bins searched), it was flagged as a 
burst. 

We found six bursts in RXTE data of 
Swift J1822.3— 1606. The burst properties are summa- 
rized in Table O In the Table are the MJDs of each 
burst, the number of counts in a 31.25-ms bin, and the 
probability that the burst would occur by chance given 
the local mean count rate. An insufficient number of 
bursts was detected to perform a detailed statistical 
analysis of the burst properties for Swift J1822.3— 1606. 
The bursts found were very narrow, typically only 
one or two 31.25-ms bins wide, and not very fluent 
compared to typic al magne tar burst s (see 'Gogu s et al 
20011: iGavrifl et al.H 2004: Scholz_&K_aspi 2011; Lin et al 



Table 5 

X-ray Bursts from Swift J1822.3-1606. 



20111) . No significant changes in the long-term flux 
decay were observed at the times of these bursts. 

Although in certain Swift observations we detected sev- 
eral bursts, these had much softer spectra than typical 
magnetar bursts and were also seen in the background 
region. Therefore, we do not believe they originate from 
Swift J1822.3— 1606. No other bursts were detected in 
any of the Swift data. 

4. DISCUSSION 



RXTE Obsid 



MJD Total counts Chance Prob"^ 



RXTE bursts 



D96048-02-01-01 


55761.53224 


15 ±4 


7.8 X 10" 


-7 


D96048-02-01-01 


55761.57082 


36 ±6 


8.6 X 10" 


33 


D96048-02-01-02 


55762.49919 


21 ± 5 


1.1 X 10" 


13 


D96048-02-03-04 


55777.91627 


12 ±3 


4.5 X 10- 


-5 


D96048-02-04-01 


55782.53122 


13 ±4 


2.4 X 10- 


-5 


D96048-02-05-01 


55789.96209 


11 ± 3 


2.2 X 10- 


-4 



Note. — ^Thc probability of the detected signal being dnc to 
noise. 

We have presented Swift, RXTE, Chan- 
dra observations following the discovery of 
Swift J1822.3-1606 during its outburst in 2011 July. 
We presented a phase-connected timing solution which 
suggests a spin-down inferred B ^ 5x 10^^ G, the second 
lowest measured for a magnetar thus far, although we 
note that timing noise may significantly contaminate 
this estimate. The flux of the magnetar was found to 
be decaying, both in total and pulsed flux, according 
to a double-exponential model. The spectrum softened 
following the outburst. Swift J1822.3— 1606 also emitted 
several short bursts during its period of outburst. We 
also a nalysed an ar c hival ROSAT observation from 
which E sposito et al.l (|2011[ ) previously reported that 
Swift J1822.3— 1606 is detected in quiescence. We note 
that the source had a similar absorption column density 
to the nearby Galactic Hii region M17. In the following 
we discuss the above results. 

4.1. Timing Behavior 

In Section 13.21 we presented a timing solution for 
Swift J1822.3— 1606 with just i' and v fitted (Solution 
1). However, this solution appears significantly contam- 
inated by timing noise, a common phenomenon in pul- 
sars. Most pulsars seem to display some unexplained 
'wandering' in their spin evolution (Hobbs et al. 2010). 
A measure of the amount of timing noise displayed by a 
pulsar is Ag and is d efined as Ag = log[(l/ 6 t^)|i>|( 10^s)'^] 
i|Arzoumanian et al.l [TOOl . iHobbs et all ()201ClD mea- 
sured a correlation between P and Ag using timing so- 
lutions for 366 rotation-powered pulsars. Magnetars are 
very noisy timers, generally having more timing noise, 
as measured by Ag , than those rotation-powered pulsars 
of sim ilar properties (jGavriil fc Kaspiir2002l : iWoods et al.l 
I2002f) . Here, for Swift J1822.3-1606, we measure Ag = 
2.8 (using i) from Solution 3) which is muc h higher than 
the va lue predicted from the correlation in IHobbs et al.l 

([2OT0h of 2. However, we caution that in general the 

i) used to calculate Ag is measured for a data span of 
10® s, whereas our data span in much shorter. The large 
value of Ag we measured may be biased by the short 
span, or by unmodelled relaxation following a hypothet- 
ical glitch that could have occured at the BAT trigger. 
Glitches are commonly seen to accompany radiative out- 
burst s from magnetars (e.g. iKaspi et al.ll2003l : iDib et afl 
12001 . 

Due to the presence of timing noise, we take the tim- 
ing and derived parameters of Solution 3 not as the 'true' 
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2-10 keV Unahsorbed Flux 



Figure 7. Spectral hardness as a function of unabsorbed flux 
for Swift J1822.3-1606. Tiie Swift WT mode observations are 
sliown as grey crosses, tlie Swift PC mode data are denoted by open 
circles, and the blue triangles respresent the Chandra observations. 



spin-inferred values, but as a 'best guess' given the data 
thus far. As such, the uncertainties in the parameters 
presented, which do not take into account the effect 
of contamination by timing noise, hkely underestimate 
the true uncertainties. Further timing observations of 
Swift J 1822.3— 1606 will help to average over the effects 
of timing noise and thus provide improved estimates of 
the spin-inferred magnetic field of the pulsar. 

The B-field measured by Solution 1 would be the sec- 
ond lowest measured for a magnetar to d ate, higher 
than only SGR 0418-1-5729 (|Rea et al.l[20Tol ). Solution 
3, although still the second lowest yet measured, gives 
a higher value of B that is close to that of magne- 
tar IE 2259-1-586 and the magnetically active rotation- 
powered pulsar PSR J1846— 0258. It is also similar 
to the quantum critical fiel d of Bqed = 4.4 x lO^'^ G 
([Thompson fc DuncanI [T996h which has been viewed in 
the past of being a lower limit on the magnetic field of 
magnetars, although SGR 0418-1-5729 has shown that it 
is not a necessary condition for magnetar-like activity. 

4.2. Flux and Spectral Evolution 

In the twisted-magnetosphere model of magnetars, the 
thermal emission is thought to originate from heating 
within the star, cause d by the decay of strong internal 
magnetic fields fiThompson et al.]l2002l ). Currents in the 
magn etosphere, which are due to twists in the m agnetic 
field ()Thompson et al.|[200l[Beloborodovl[2009l) . scatter 
the thermal surface photons to higher energies. In addi- 
tion to scattering, the currents provide a source of sur- 
face heating in the form of a return current. The flux 
increase that accompanies a magnetar outburst is the- 
orized to be due a rapid heating which could originate 
from magnetospheric, internal, or crustal reconfiguration 
of the neutron star. This release may result in a signif- 
icant increase in the surface temperature, in the return- 
current heating, and in the twisting of the magnetic field. 
Thus, an increase in flux due to an internal heat release 
should result in an increase of the hardness of the emis- 
sion. This hardness-flux correlation is in agreeme nt with 
observations of severa l magnetars ([G otthclf & HalperE 
2007t iTam et al.|[2008t iZhu et al.|[2008t iScholz fc Kaspi 
2011 ). 



IScholz fc Kaspil (|2011f) explored the hardness-flux cor- 
relation for magnetar outbursts by comparing the rela- 
tion between fractional increase in 4-10/2-4 keV hard- 
ness ratio and fractional increase in 2-10 keV unabsorbed 
flux for six different outbursts in four different magne- 
tars. We present a similar plot here in Figure [7] for 
Swift J1822.3— 1606. Here, however, the hardnesses and 
fluxes are absolute quantities and not fraction a l incre ases 
over quiescent values as in IScholz fc Kaspil ()201lD . as 
there is no quiescent observation of Swift J1822.3— 1606 
with the appropriate spectral coverage. Figure [7| shows 
that Swift J1822.3— 1606 softens as the flux decreases fol- 
lowing the outburst and so is in broad agreement with 
the hardness-flux correlation observed in other magnetar 
outbursts. This spectral softening with flux decline is 
clear also in Figure El where kT declines and the power- 
law index F increases as the flux drops. 

4.2.1. Magnetars in quiescence 

The quiescent flux of Swift J1822.3— 1606 mea- 
sured by ROSAT in 1993 is about 3 orders of 
magnitude lower than the peak flux measured fol- 
lowing the outburst. Such large flux variations 
have been observed in several other magnetars (e.g. 
IE 1547-5408, XTE J1810-197, AX J 1845-0258; 
Ibrahim et all 120041: [Gotthelf et al.' '2 0041: iTam et al.l 



20061 : IScholz fc Kaspil 120111 : iBernardi niet alJ 120111). 
Othe r magnetars , such as IE 1841—045 fZhu & Kasgi 
2OT0t 'Lin et al] l20Tl[ ). 4U 0142-^61 (iGonzalez et al 
and IRXS J170849. 0-400910 (Iden Hartog et al 



2008) have not exhibited large flux variations, but are 
much brighter in quiescence than are the magnetars with 
large outbur s ts. T he cause of this difference is unclear. 
iPons fc Real (|2012D suggest that there IS a maximum 
luminosity that can be reached by a magnetar during 
an outburst due to neutrino cooling dominating at high 
crust temperatures. This helps to explain the differences 
in outburst magnitudes, but does not address the wide 
range of quiescent luminosities. 

Case in point, the magnetar IE 2259-1-586 has spin 
properties that are likely quite similar to those of 
Swift J1822.3— 1606 but has a much higher quiescent lu- 
minosity. The magnetic field me asured from spin-dow n 
for IE 2259-^586 is 5.9 x 10" G (jGavrifl fc Kaspil[200l . 
close to B = 5.1 X lO^^ G for Swift J1822.3-1606 as es- 
timated by our Solution 3. IE 2259-1-586 also went into 
a period of outburst on 200 2 June 18 where t he flux 
increased by a factor of ^ 20 (jWoods et al.ll2003 l. How- 
ever, in quiescence, IE 2259-1-586 is much brighter than 
Swift J1822.3— 1606 with a quiescent 2-10 keV luminos- 
ity of 2 X 10^4 erg s-i (jZhu et all [20081 ) compared to 
< lO-'^^erg s-i for Swift J1822.3-1606. 

One possibility is that the 'true' magnetic fields of the 
more luminous magnetars are higher than those of the 
fainter magnetars. The spin-down of the neutron star 
is only sensitive to the dipole component of the mag- 
netic field. If the magnetic field had significant com- 
ponents in higher multipole s or a toroidal comp onent 
ijThompson fc DuncanI [l996l : iPons fc Perilal l20Tll ). the 
true magnetic field could be higher. 

Another possibility is that neutrino cooling in the 
core is setting a long-term luminosity limit, and that 
the neutrino cooling properties of the stars are differ- 
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Figure 8. A model of the thermal relaxation of the neutron- 
star crust that approximately reproduces the observed 1-10 keV 
luminosity decay assuming a distance of 1.6 kpc. In the model, we 
deposit an energy 3 X 10*^ ergs in the outer crust between densities 
of 2 X 10^-3 X 10^" g cm~^. The subsequent cooling of the crust 
gives a good match to the observed light curve. 

ent, e.g. due to different masses. For example, con- 
sider first the case where the neutrino emission in the 
core is due to the modified URCA proces s, with an 
emiss ivity ~ 10^° erg cm~^ s~^ Tg ( Yakovlev et al.l 
|2003() . If we take the magnetic- field decay time to 
be T = 10'* yrs, then the luminosity from magnetic 
field decay is roughly Lb = {B"^ /STr){47rR^ /3){1/t) = 
10^"* erg s~^ for B = 10^"* G. Balancing this with the 
neutrino losses L^, — {AttR^ /3)e,^, we find a core temper- 
ature Tc — 2.5 X 10* K or, using the core temper atu re- 
luminosity relation from iPotekhin fc YakovlevI (|2001h , a 
luminosity L « 4 x 10^'^ erg s~^. On the other hand, 
if the neutrino emission is by the direct URCA process , 
with - 10^6 erg cm'^ s'^ r| (| Yakovlev et al.l[200l . 
we find a core temperature Tc = 1.5 x 10^ K, correspond- 
ing to a surface luminosity of w 2 x 10^^ erg s~^. This 
shows that we might reasonably expect a factor of > 200 
in luminosity between different stars if one has slow neu- 
trino emission in the core, and the other fast, for example 
if the mass of one of the stars is large enough for direct 
URCA reactions to occur in the core. Even in the case 
where external currents dominate the quiescent luminos- 
ity, thermal emission from the neutron star provides a 
baseline luminosity, so that the low quiescent luminosity 
of Swift J1822.3— 1606 suggests a low core temperature 
which implies either a low heating rate or efficient neu- 
trino emission. 

4.2.2. The observed luminosity decay of Swift ,71822.3-1606 

We find that the observed luminosity decay is well re- 
produced by models of thermal relaxation of the neutron- 
star crust following the outburst. An example is shown 
in Figure O which shows the cooling of the crust after 
an injection of « 3 x 10**^ ergs of energy at low density 
« 10 ° g cm~'^ in the outer crust at the start of the out- 
burst. We follow the evolution of the crust temperature 
profile by integrating the thermal dif fusion equation. The 
calculation and microphysics follow [Brown &: Cummin^ 
(|2(309,;) who studied transiently accreting neutron stars, 
but with the effects of stron g magnetic fields on the 
thermal conductivity included (jPotekhin et al.|[T999l) and 
for the outer boundary condition using the Tee-T-mt re- 
lation appropriate for a magnetized envelope following 
iPotekhin fc Yakovlev, (,2001.) . The calculation follows the 



radial structure only; we assume that the magnetic field 
geometry is a dipole and take appropriate spherical av- 
erages to account for the variation in the rmal conduc- 
tivity across the star (|Potekhin et al.iri999f ). We assume 
i? = 6 X 10^^ C, similar to the value inferred from the 
spin down, a 1.6 Mq, i? = 11.2 km neutron star, and take 
an impurity parameter for the inner crust of Qimp = 10 
()Jonesll20"M ). We set the neutron-star core temperature 
to 2 X 10^ K, which is needed to obtain a quiescent lu- 
minosity < 10^^ erg s"""^. 

With the neutron-star parameters fixed, we then vary 
the location and strength of the heating and find that we 
obtain good agreement with the observed light curve for 
times < 100 days, if the initial heating event is located at 
low densities < 3 x 10^*^ g cm~'^. This conclusion comes 
from matching the observed timescale of the decay, and 
is not very sensitive to the choice of neutron-star param- 
eters. For example, changing the neutron-star gravity 
changes the crust thickness and therefore cooling time, 
giving an inferred maximum density pmax g~^. This 
means that the inferred density can change by a factor of 
a few but cannot be moved into the neutron drip region, 
for instance. That only a shallow part of the outer crust 
is heated is an interesting constraint on models of crust 
heating in a magnetar outburst. We find that it is dif- 
ficult to match the observed light curve at times > 200 
days, but the late time behaviour of the light curve is 
sensitive to a number of physics inputs associated with 
the inner crust, including the thermal conductivity and 
superfluid parameters, as well as modification due to the 
angular distribution of the heating over the surface of the 
star. We will investigate the late-time behaviour in more 
detail in future work. Figure |8] suggests that the source 
could undergo significant further cooling in the coming 
years. 

4.3. Distance Estimate and Possible Association 

As shown in the ROSAT image (Figure [T]), the 
Galactic Hii region M17 is located ~20' southwest of 
Swift J1822.3-1606 It has a distance of 1.6 ± 0.3 kpc 
(|Nielbock et al.|[200ll ) and an absorption column density 
A^H = 4 ± 1 X 10^1 cm -2 (jTownslev et al.ll200l which 
is consistent with our best-fit value of 4.53 ± 0.08 x 
10^1 cm -2. This suggests that Swift J1822. 3-1606 could 
have a comparable distance to that of MITFI. If so, then 
Swift J1822.3— 1606 would be one of the closest magne- 
tars detected thus far. 

The above argument does not necessitate a direct as- 
sociation between M17 and Swift J1822.3— 1606. How- 
ever, if Swift J1822.3— 1606 is associated with M17, then 
its angular separation of 26' from the cluster center, 
where the X-ray emission peaks in the ROSAT image, 
implies a physical distance of 12 pc. For a pulsar age of 
10^ yr, this requires a space velocity of only r-^100kms~^ 
(corresponding to a proper motion of 0'.'016yr~*). This 
would make a direct proper motion measurement diffi- 
cult. From timing, the characteristic age appears to be 
larger than 10^ yr which would further reduce the im- 
plied proper motion. On the other hand, characteristic 

While there are two molecular clouds surrounding M17 
llWilson et al.l [20031 ). they are confined to the north and west, 
such that they should not contribute to the TVh of either M17 
or Swift J1822. 3-1606. 
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ages can be large overestimates of the true age. However, 
even if the true age were as low as 10^ yr, the proper mo- 
tion would be difficult to measure even with Chandra. 
Additionally, if the magnetar was born near an edge of 
the cluster, the angular separation from its birthplace 
could be larger or smaller by up to ~ 10'. 

5. CONCLUSIONS 

We have presented an analysis of the post- 
outburst radiative evolution and timing behavior of 
Swift J1822.3— 1606, following its discovery on 2011 July 
14. Following a timing analysis for the source post- 
outburst, we estimate the surface dipolar component of 
the B-field to be 5 x 10^^ G, slightly higher than that 
inferred in Paper I. However, as this measurement is con- 
taminated by timing noise, the true value of the magnetic 
field could be well outside of the uncertainties quoted in 
Table [3 Futhcr monitoring of Swift J1822. 3-1606 as 
it fades following the outburst will allow us to better 
account for the timing noise and measure more robust 
timing parameters. 

The quiescent flux of Swift J1822.3— 1606 measured 
using a 1993 ROSAT observation of Ml 7 was found 
to be roughly three orders of magnitude lower than 
the peak flux during the outburst. The flux evolu- 
tion following the outburst was well characterized by 
a double-exponential decay. By applying a crustal 
cooling model to the flux decay, we found that the 
energy deposition likely occured in the outer crust at 
a density of ^ 10^°gcm~'^. The spectral properties of 
Swift J1822.3— 1606 were observed to soften following 
the outburst, with the power-law index increasing and 
the temperature of the blackbody decreasing. Indeed, 
a hardness-flux correlation, similar to what is observed 
in other magnetars, was clearly observed. Based on the 
similarity in Nh to that of the Hii region M17, we argue 
for a source distance of 1.6 ± 0.3 kpc, one of the closest 
distances yet inferred for a magnetar. 
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